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We study the role that phase disorder plays in the mode discrimination characteristics of a Talbot res- 
onator. Two cases considered here correspond to a six-core and a twelve-core photonic crystal fiber. 



I. INTRODUCTION 

Fiber laser arrays are currently considered as one of the most viable options for state of the art high power 
lasers [l| . Key to their effectiveness is to coherently combine the output of as large a number of amplifiers as 
possible. Coherent beam combining can be implemented either passively 043 , or actively @, with the first 
one being preferred given that in principle it would be easier to implement. Passive schemes are based on 
linear mode coupling of optical elements and with the advancement of technology in multicorc fiber arrays, 
many of the recent schemes make use of such arrays. Regardless of the beam-combining architecture in all 
instances, we define E output = J2i J2j ^ijEj as the output field that results from combining electric fields 
Ej emerging from individual and ideally identical amplifiers. Here Sy defines the passive coupling scheme. 
The objective of any combining scheme is effectively one of efficient mode selectivity or mode discrimination. 
This is the case if for example there is an eigenmode of Sy whose (largest) eigenvalue is much bigger than 
the rest. Furthermore maximum coherence is achieved if the dominant mode is such that the phases of 
j^{out) _ ^ SijEj are identical. Having said this, it is typical that in most if not all instances, as the 
number of elements N in the arrays increase, the coherence decreases [Io[. This is due to the fact that 
most of perturbations appearing during the propagation of the fields in the (ideally identical) individual 
amplifiers affect the path length, thus in turn the phase. This combined with the fact that as N increases 
mode discrimination decreases, one expects coherence to diminish. For each passive linear scheme, it is 
important to characterize its robustness to phase disorder; in particular it would be useful to determine how 
the coherence changes when there is disorder in the input phases and how it scales with N. 

Given the fast technological advances in multicore fibers, they present perhaps the most suitable archi- 
tecture for passive beam combining and in this paper we discuss such multicore arrays in a Talbot resonator 
fill ] as proposed in [l2|. The natural modes in this resonator are due to self- imaging, which according to 
the Talbot effect happens to images formed by periodic objects, such as a periodic diffraction grating. The 
combination of periodicity and Frcsncl diffraction is such that the image (mode) is reproduced every multiple 
of the Talbot length zt — , where a is the period of the grating and A is the light wavelength. 

The next section presents the main results highlighting the role of phase disorder on mode selectivity. 
While our work deals with one particular beam combining scheme, we believe this behavior is generic of all 
passive, linear schemes. This is why in section HTT1 we briefly discuss nonlinear schemes and their suitability 
to achieve high coherence for large arrays. 
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II. TALBOT RESONATOR 



In this section, we present the results of numerical simulation about the role of phase disorder on beam 
combining. We limit our study to the scheme based on a Talbot resonator for a photonic crystal fiber 
geometry fl2| . We will consider two periodic like profiles: a hexagonal one consisting of six cores, equally 
spaced, one at the center of the hexagon. The second case is a rectangular array of 12 cores. In either case 
the generated modes will take the form 



E m (x,y,0) = ]Texp(-a[(* - x 3 f + (y - y 3 ) 2 ])e^ m> (1) 

with n = 6 or n = 12. 

The figure of merit to be studied is the correlation function, 



/ s = \J E m(x,y, z = 0)E m (x,y,z)dxdy\ 
1{Z) f\E m (x,y,z = 0)\2dxdy ' U 

where ideally you would like to work in a regime for which 7 is not only high for the in-phase mode, but 
separates well from other modes. These two properties are, for the Talbot resonator what determines if 
mode selectivity is achieved. 

Figures 11121 show the dependence of correlation function 7(2;) on distance from the output of waveguide in 
a 6-core and 12-core PCFs, respectively. The distance is in units of the lattice constant A, where A/A = 4.51. 

For each supermode to, we consider the deterministic phase factor <j^™^ = 2m(j — l)-7r/6 associated with core 
j, where j = 1, • ■ ■ ,6 for 6-core PCF and j = 1, • ■ ■ ,12 for 12-core PCF. As shown in both figures \T\2\ it is 
the case that in the absence of disorder both the N = 6 cores in a hexagonal lattice and the N = 12 cores in 
a rectangular lattice, mode selectivity is achieved. We should point out that figure [T] reproduces the figure 
9 in [Hi. 

The rest of the figures represent averaged values (7) computed over many realizations of different strengths 
of the random phase perturbations. We assume in all instances the phase inputs are iid. Figures 02 |4j [5] 
correspond to different realizations of the random phase factor in a 6-corc PCF with same parameters as 
figure [1] for increasing disorder strengths. That is, in our simulations the phase includes a random variable £ 
uniformly distributed so that (f)j — 2m(j — l)7r/6 + £ and with strengths < £ < tt/2 (figure|3|), < £ < tt 
(figure EQ) and < £ < 2tt (figure [5]). We clearly observe that while the contrast (separation of different 
curves) between the in-phase (m = 0) mode and the others does not diminish by much for noise levels 
< £ < 7r/2 (figure 3) it docs so in a significant way once the noise strength increases (figure 4). That 
is, even if one were able to control path lengths so that the noise is only half a period < £ < tt this 
scheme becomes ineffective. Notice that if < £ < 2ir (figure [S]) there is no mode discrimination at all. We 
should notice that for the operational wavelengths in fiber lasers, small perturbations (e.g., nonlinear phase, 
fiber lengths, thermal effects, gain induced perturbations to name some ) are likely to modify the phase in 
multiples of 2ir. Figure |5] shows the average of 7(2) over different realizations of random phase factor in a 

12-core PCF, where (j) ( - n) = 2m(j - 1)tt/3 + £ for j = 1, ■ ■ ■ , 12, m = 0, 1, 2, and < £ < 2tt. The similarity 
of the outcome with the 6-core case suggests that noise will "erase" the contrast irrespective of the number 
of fibers or the geometry. 
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FIG. 1: Dependence of correlation function 7(2) on distance from the output of waveguide in a 6-core PCF. The 
phase factor <?^ m ' is considered deterministic. The distance z is in units of lattice constant F, where A/A = 4.51. 
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FIG. 2: Dependence of correlation function 7(2) on distance from the output of waveguide in a 12-core PCF. The 
phase factor <j>^ m ' > is considered deterministic. The distance z is in units of lattice constant F, where A/A = 4.51. 



III. NONLINEAR COUPLING SCHEMES 



The previous section highlights the sensitivity to phase disorder and the loss of coherent beam combining 
efficiency. Furthermore, it is generally the case (with or without disorder) that as the number of elements in 
a fiber laser array increases, the coherence diminishes for linearly based coupling schemes [Ioj |. The role of 
phase disorder in linearly based coupling schemes we believe will be similar to what we just showed forthe 
Talbot resonator discussed above. Altogether, new paradigms need to be sought. Here we suggest that 
having a nonlinear coupling scheme should be considered as a viable alternative. As we know from many 
studies not only in laser dynamics but in other areas of science, synchronization of weakly coupled oscillators 
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FIG. 3: Average of 7(2) for N = 6. Random phase factor is considered as </>' m ' ) = 2m(j — + £ for < £ < tv/2. 
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FIG. 4: Average of 7(2) for N = 6. Random phase factor is considered as <K = 2m(j — l)-7r/6 4- £ for < £ < 7r 



is universal. Whether one describes laser modes, social networks, fireflies or chemical reactions to name some, 
a phase transition where incoherent oscillations lock into an in-phase state is present at a critical coupling 
strength. For this to happen, the system requires to be nonlinear, but the details need not be identical. This 
behavior is best represented in the work by Kuramoto on a "globally" weakly coupled nonlinear chain of 
oscillators, 



JV 



^ = + If Yl sin (^ _ & * = !> 2 > •••> N ' 



(3) 



i=i 



where 9i (£) is the phase of the i — th oscillator and w.; is its natural frequency. In this model, disorder comes 
from the fact that the frequencies are random with a Lorenzian probability density 
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FIG. 5: Average of 7(2:) for N — 6. Random phase factor is considered as 4>j = 2m(j — l)7r/6 + £ for < £ < 2ir. 




50 100 150 200 

distance (A) 



FIG. 6: Average of 7(2) for N = 12. Random phase factor is considered as <f>j = 2m(j — 1)tt/3 + £ for < £ < 2n. 



= 1x2 : / \2r ( 4 ) 



Kuramoto [13j found a phase transition in that the order parameter 



1 N 

r(*) = i^E e ^ (i) i ( 5 ) 
3=1 



takes the asymptotic (N — > 00 and i — > 00) value 



G 



r = 0, K < 26 (6) 

which translates into enhanced locking as the coupling strength goes past a critical value. 

Perhaps more relevant case to coherent beam combining which can still be viewed as some form of phase 
transition is that of spatial and spatio-temporal localization in fiber arrays. Here the nonlinear index of 
refraction together with linear nearest neighbor coupling is modeled by the discrete nonlinear Scrodinger 
equation (DNLSE), 



.da r , 



i- 



(Aa) nm + n 2 \a n m\ a nm , (8) 



where the term Aa represents a general linear coupling scheme. Specific to cw lasers, the possibility of 
reaching an output that is spatially localized into a few cores would enhance coherence if it is measured by 



As theory 14 1 and experiment [15| have shown, at sufficiently high powers (the equivalent of having high 
enough K in the oscillator model) , an initially broad distribution of power reaches in propagation a localized 
mode. We are currently exploring the role that the disorder of incident phases plays in the dynamics of 
spatial and spatio-temporal localization. It is also the case that non Hamiltonian perturbations to coupled 
systems like the DNLSE may produce an attractor in the form of a highly coherent state. As a proof of 
principle, recent work fliij on large systems of nonlinearly coupled oscillators shows one can find asymptotic 
dynamics towards an attracting highly coherent state. Figure (J7J shows a dynamics for an array of 17 
oscillators reaching a state of maximum coherence. 

Finally, nonlinear coupling schemes can perform in a very robust manner even in the presence of losses. In 
fact, small losses (that is departing from a conservative system where stable equilibrium states are centers) 
should facilitate synchronization |17| . 



IV. CONCLUSIONS 



Phase disorder is always present in the amplification and propagation of light in optical fibers. This work 
presents the deterrent effect it has on the efficient beam combining in a Talbot resonator. Two cases of 6 and 
12 elements are considered as a way to explore the role if any, N the number of elements plays on the efficiency 
of coherent coupling. We found for these two cases that mode selectivity characterized by the separation 
of the parameter 7 for different modes lost if the noise in the phases is of strength equal to 27r. While 
the results are limited to one particular linear scheme, we believe all linear schemes will show an eventual 
decrease in efficiency not only by increasing the number of elements as previously stated [Ioj |. but also in 
terms of expected phase disorder generated during propagation in the individual fiber amplifiers. The last 
section provides a brief discussion of a new paradigm whereby considering nonlinear coupling mechanisms, 
phase transitions towards high coherence are at least theoretically feasible in the important N — > 00 case. 
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FIG. 7: Time evolution of E = J^i , for an initial condition with E — 0.528646241 for an array of 17 nonlinear 
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coupled oscillators model [16| 
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